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ABSTRACT 

Condensate clouds are a salient feature of L dwarf atmospheres, but have been assumed to play 
little role in shaping the spectra of the coldest T-type brown dwarfs. Here we report evidence of 
condensate opacity in the near-infrared spectrum of the brown dwarf candidate Ross 458C, obtained 
with the Folded-Port Infrared Echellette (FIRE) spectrograph at the Magellan Telescopes. These data 
verify the low-temperature nature of this source, indicating a T8 spectral classification, log 10 L^oi/Lq 
= — 5.62±0.03, T e ff = 650±25 K, and a mass at or below the deuterium burning limit. The data 
also reveal enhanced emission at i^-band associated with youth (low surface gravity) and supersolar 
metallicity, reflecting the properties of the Ross 458 system (age = 150-800 Myr, [Fc/H] = +0.2 to 
+0.3). We present fits of FIRE data for Ross 458C, the T9 dwarf ULAS J133553.45+113005.2, and the 
blue T7.5 dwarf SDSS J141624.08+134826.7B, to cloudless and cloudy spectral models from Saumon 
& Marley. For Ross 458C we confirm a low surface gravity and supersolar metallicity, while the 
temperature differs depending on the presence (635^5 K) or absence (7601™ K) of cloud extinction. 
ULAS J1335+1130 and SDSS J1416+1348B have similar temperatures (595±H K), but distinct surface 
gravities (log<? = 4.0-4.5 cgs versus 5.0-5.5 cgs) and metallicities ([M/H] « +0.2 versus -0.2). In all 
three cases, cloudy models provide better fits to the spectral data, significantly so for Ross 458C. 
These results indicate that clouds are an important opacity source in the spectra of young cold T 
dwarfs, and should be considered when characterizing the spectra of planetary-mass objects in young 
clusters and directly-imaged exoplanets. The characteristics of Ross 458C suggest it could itself be 
regarded as a planet, albeit one whose cosmogony does not conform with current planet formation 
theories. 

Subject headings: stars: fundamental parameters stars: individual (Ross 458C 

(ULAS J130041. 72+122114.7), ULAS J133553.45+113005.2, 

SDSS J141624.08+134826.7B (ULAS J141623.94+134836.3)) stars: low 
mass, brown dwarfs — stars: planetary systems 



1. INTRODUCTION 

Mineral condensate clouds are a unique and promi- 
nent constituent of the atmospheres of the coldest 
brow n dwarfs, the L dwarfs and T dwarfs (Kirkpatrick 
120051 and references therein). For the L dwarfs, con- 
densates are a primary contributor to shaping spec- 
tral energy distributions (SEDs). Their formation de- 
pletes the atmosphere of TiO and VO gases, trans- 
forming optical SEDs and defining t he L dwarf spec- 
tral s equence dJones fc Tsuiil 119971 : iKirkpatrick et al.l 
Il999t lLoddersj I2002D . Condensate opacity contributes 
to muted infrared H2O bands and the characteris- 
tic red near-infrared colors of L dwarfs (J — K = 
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1.5- 2. 5; IChabrier et all 1200(1 IKnann etail l200l . Ev- 
idence of silicate grain absorption has al so been found 
in the mid-infrared spectra of L dwarfs (iCushing et all 
[20061: iHelimg et al.1 [20061: ILooper et all I2008D . For the 
T dwarfs, condensates appear to play a minor role, 
as evidenced by strong molecular gas absorption fea- 
tures and typically blue near-infrared co l ors (J — K = 
-O.5.. 0.5: IBurgasser et all l2002al [200611 iGeballe et all 
[2001 iCushing et all 12005ft . This shift is believed to 
be the result of clouds sinkin g bel ow the visible pho- 
tosphere (lAckerman fc Marlevll200TI: iCooper et alll2003 ; 
i Tsuiil 120051) or breaking apart (jBurgasser et al . 2002b; 
iMarley et all Isubmittedft , although the details of this 
transition remain poorly understood (Burrow s et all 
2006; IBurgasserl[2007l: ISaumon fc Marlevll2008ft .~ 

Even without considering condensates, modeling the 
highly structured SEDs of T dwarfs remains a chal- 
lenge. The abundances and opacities of pro minent 
gas s pecies - H 2 (collision induced absorption; iLinskyl 
1969), H 2 0, CH4, NH 3 and neutral alkalis (p ressure- 
broadened wings ; iBurrows fc Volobuvevl |2003) - are 
highly sensitive to photospheric gas conditions. Tem- 
perature and pressure variations, elemental abundances, 
and vertical mixin g: all have observable effects on T 
dwarf spectra (e.g.. | Knapp et all 120041: IBurgasser et all 
I2006al: ISaumon et all 120071: iStephens et all 12009ft . Li 
principle, one could disentangle the contributions of 
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these atmospheric properties on observed spectra and 
derive bulk characteristics — ma ss, age and meta l licity— 
for individua l sou rces (e.g., Burgass er et al.l [2006a; 
iLeggett et al.1 l2009f >. However, a variety of effects 
drive systematic deviations between model and ob- 
served spectra. These include incomplete opacities 
for key abso rbers such as CH4, NH? and K I (e.g., 
Allard et al.l l2003t iCushing et all 1200a iFreedman etaJ] 
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H); non-e quilibrium chemistry arisin g from vertical 
mixing (e.g., iSaumon et al.l [20061 l2007f): metallicity ef- 
fects on chemical pathways (e.g-JLodders fc Fegleyll2002l ; 
IBurgasser etaHl200a IVisscher et alJl2010D : and the in- 
fluence of conde nsate clouds in regions of minimum 
gas opacity (e.g.jTsuii et"aLlll999t lAckerman fc Marlevl 
2001 1 ; IBurrows et al.ll2002F These processes lead to per- 
sistent uncertainties in inferred atmos pheric and physical 
parameters for individual sou rces (cf.. lBurningham et al.1 
2D1H: IBurgasser efall 1201 Obi) . 

Benchmark T dwarfs, with independent determina- 
tions of age, mass and metallicity, are imp ortant calibra- 
tors f or brown dwarf atmospheric studies dPinfield et al.l 
2006). Such benchmarks include widely-separated 
T dwarf companions to nearby stars, several exam- 
ples of wh ich have been identi f ied over the past 15 
years (e.g.. iNakaiima et aT] Il995t IBurgasser et al.l [2~000t 
iLuhman et al.ll2007t iBurningham et all 120091). Mos t re- 
cently. IGoldman et al.l (|2010l) and iSchold (|20"T0al) re- 
ported the discovery of a faint (J = 16.67±0.01), 
co-moving candidate T dwarf companion to the Ross 
458 system in the UKIRT In f rared Deep Sky Sur- 
vey (UKIDSS; lLawrence etHI 12051 . The source, 
ULAS J130041. 72+122114.7 (hereafter Ross 458C), is 
separated by 102" from the M0.5 Ve + M7 Ve Ross 
458AB pair, implying a projected separ ation of over 
1100 AU at a distance of 11.7±0.2 pc (jvan Leeuwenl 
120071 ). The late- type nature of Ross 458C was in- 
ferred from its optical and near-infrared colors (z — J = 
3.55±0.19; J - K = -0.21±0.06), faint absolute magni- 
tude (Mj = 16.40±0.04) and large fl ux contrast in filter s 
sampling the 1.6 /jm CH4 band (cf.. lTinnev et al1[ 2005). 
Fro m these photom e tric m easurements, iSchold (|2010af ) 
and IGoldman et al.1 (|2010l ) estimated spectral types of 
T7±l and T8-9 for Ross 458C, respectively. However, 
no spectrum was reported in either study. 

If Ross 458C is a brown dwarf, it is a potentially 
interesting benchmark for several reasons. Kinemat- 
ics of the R oss 458 system coincide with the Hyades 
supercluster (jEggenl [l960t iMontes et al.1 120011 ). an ex- 
tended population of relatively young and metal-rich 
star s that includes the 62 5±50 Myr Hyades open clus- 
ter (|Perrvman et al.lll998l ). The prima ry exhibits con- 
spicu ous activity in Ha and Ca II lines ([B rowning et al. 
20101). and is an X-ray sou rce (jKiraga fc Steoied 120071 
Lopez-Santiago et al.l [2009). Ross 458C could there- 



fore be a young and/or metal- rich brown dwarf bench- 
mark, useful for constrai ning both evolutionary and 
atmospheric model s (cf. iZapatero Osorio et al.l 120041 : 
IDupuv et alll2009al ). 

In this article, we present the first spectroscopic obser- 
vations of Ross 458C, obtained in the near-infrared with 
the ne wly-commissioned F olded Port Infrared Echellette 
(FIRE: lSimcoe et al.ll2008h . These data confirm the late- 
type T dwarf nature of the source, and indicate a low 
surface gravity and supersolar metallicity consistent with 



characteristics of the primary. Furthermore, model fits 
indicate that cloud opacity plays a prominent role in 
shaping its near-infrared spectrum. In Section 2 we con- 
strain the age and metallicity of the Ross 458 system 
based on the photometric, spectral, kinematic and rota- 
tional properties of its M0. 5 Ve primary. In Section 3 we 
describe our FIRE observations of Ross 458C and two 
other late-type T dwarfs, from which we derive spectral 
types and basic physical characteristics. In Section 4 we 
fit the spectral data for these three s ources to cloudless 
and cl oudy atmospheric models from ISaumon k Mar ley 
( 2008, hereafter SM08), and compare the inferred atmo- 
spheric and physical parameters to prior studies. In Sec- 
tion 5 we discuss the role of condensates in shaping the 
near-infrared SEDs of young, cold T dwarfs, and their 
influence on T dwarf temperature scales. In Section 6 
we summarize our results. 

2. CHARACTERIZING THE ROSS 458 SYSTEM 

2.1. Examining Association with the Hyades Open 
Cluster 

Under the assumption of coevality, some physical prop- 
erties of Ross 458C can be inferred from its stellar com- 
panions, specifically age and metallicity. As such, it is 
necessary to examine the association of this system with 
the Hyades supercluster, a kinematic stream composed 
of several dozen nearby stars with similar space motions 
(Eggei3 H965t IMontes et aH 12551 . While co-moving as- 
sociations typically have similar ages and compositions, 
the ages of the Hyades supercluster stars span the star- 
formation history of the d isk (0.4-2.0 Gvn lFamaev et aT] 
120051 : lAntoia et al.l 120081 ) . so stream membership is not 
a particularly useful age indicator. Within the stream 
is the eponymous Hyades open cluster, the nearest open 
cluster to the Sun, whic hi constitutes a coeval g roup with 
an age of 625±50 Myr ([Perrvman et alJll998fT . Because 
of its proximity, members of the Hyades cluster halo are 
distributed broadly across the sky. It is therefore not un- 
reasonable for Ross 458 to be associated with this struc- 
ture. 

To examine this possibility, we calculated the systemic 
space velocity of the Ross 458 system. We used the long- 
baseline proper motion from Tycho-2 {jj, a = — 640.1±1.5 
mas yr" 1 ; (j, s = -25.1±1 .4 mas vr" 1 : iHpg et al.l l2000h. 
the radial velocity from iNidever et al.l (|2002l . V r = - 
11.2±0.1 km s _1 ), and the revised H ipparcos parallax 
(n = 85.5±1.5 mas: Ivan Lee uwen 2007). From these we 
find heliocentric (U, V, W) = (30.6, -18.8, -9.9)±(0.5, 0.3, 
0.2) km s _1 . The veloc ity of the Hyades clust er in the 
same coordinate system (|de Bruiine et al.|[2001l) is (U, V, 
W) = (42.25, -19.06, -1.45)±(0.14, 0.13, 0.21) km s" 1 , a 
diffe rence of (11.7, -0 3, 8.5 )±(0.5, 0.3, 0.3) km s" 1 (see 
also IPerrvman et alj 119981 ) . This is much larger than 
the ID velocit y dispersion o f Hyades cluster members, 
-0.35 km s" 1 |MadseiJl20"5a) . amounting to a ~20a dis- 
crepancy. Hence, Ross 458 is not currently kinematically 
associated with the Hyades open cluster. 

Has Ross 458 ever been associated with the Hyades 

9 We adopted a heliocentric, right-handed rectangular coordinate 
system with U velocity locally pointing away from the Galactic 
center, V velocity pointing in the direction of Galactic rotation (I, 
b = 90°, 0°) and W velocity pointing towards the north Galactic 
pole (6 = +90°). 
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cluster? Its current separation from the cluster center is 
53 pc. Using the epicyclic orbit approximation to model 
the past orbits of Ros s 458 and the Hyades (a dopting 
the Oort parameters of Fe ast fc Whitelocklll997l and Lo- 
cal Standard of Rest from iDehnen fe Binnevl[l998l ). we 
find that Ross 458 was never closer to the Hyades than 
~50 pc during the past 100 Myr, or about 5 tidal radii. 
Hence, there is no evidence of past or present physical as- 
sociation between Ross 458 and the Hyades open cluster, 
and association with the Hyades supercluster provides no 
meaningful constraint on the age of this system. 

2.2. The Age of the Ross 458 System 

Despite the wide range of ages associated with the 
Hyades kinematic stream there is considerable evidence 
that the Ross 458 sy stem itself i s qui te young. The 
M0.5 Ve primar£3 ( Reid et all H995) is an "ultra- 
fast r otator" with a perio d between 1.54 and 2.8 9 
days (jPizzolato et all 120031 : iKiraga fc Stepienl I2007T ). 
Ross 458A had the highest measured projected rota- 
tional velocity (wsini = 9.7±0.5 km s" 1 ) among 69 
nearby M0-M3 dwarfs in the California Planet Search 
program, in which only ~3% of early-M dwarfs had any 
detec ta ble rotation ( vs'mi > 3 km s _1 ; iBrowning et al.l 
120101 ). iBeuzit et al.l (|2004h have argued for an age less 
than 1 Gyr on the basis of this rotation. 

In addition, Ross 458A exhibits strong magnetic ac- 
tivity. The presence of Ha emissio n indicates an up - 
per activity age limit of 400-800 Myr ([West et alll2008l ). 
This limit is rough ly consistent with the analysis of 
iTakeda et all (|2007l ). who deduce a 68% upper limit 
of 440 Myr based on theoretical modeling of the pri- 
mary's high-resolution spectrum (however, this does not 
appear to be a "well-defined" age according to that 
study). Ross 458 is also an X-ray source, detected in 
three Rontgen Satellite (ROSAT) Position Sensitive Pro- 
portion al Counter (PSPC ) pointings and in the All-Sky 
Survey (jVoges et al.ll2000| ). Combining these four obser- 
vations yields a mean ROSAT PSPC flux in the 0.1-2.4 
keV band of 4.35 ct s _1 , and mean hardness ratios of HR1 
= -0.13 and HR2 = -0.09. The mean X -ray flux, cou- 
pled w ith the energy conversion factors of iFleming et all 
(1995) and the Hipparcos parallax, implies a soft X-ray 
luminosity of Lx = 5 x 10 29 erg s _1 and an X-ray to bolo- 
metric luminosity ratio of log Lx /Li> i = -2.7. Such a 
high ratio defines Ross 458 as a "saturated" X-ray emit- 
ter, as expected for stars of similar color and rotation 
periods of <4 days (|Pizzolato et al.1 [20031 ) . The combi- 
nation of color and rotation period for Ross 458A also 
classifies i t as a C-sequen ce ( "convective" ) rotator in the 
scheme of Barnes, (200 J). The characteristic spin-down 
timescale for such stars is ~300 Myr. From a cross- 
comparison of age-dated c luster and field populations of 
low-mass, solar-type stars. [Barnes! ([20031 ) concludes that 
the fraction of C-sequence ultrafast rotators reduces to 
zero by age ^800 Myr. Hence, this age can be consid- 
ered a robust upper limit for the Ross 458 system that 
is consistent with activity estimates. 

10 The M dwarf components of Ross 458 are separated by 0'.' 5 
IIBeuzit ct al. 2004), implying that reported optical spectroscopy 
of this system probably includes both components. However, 
Ross 458B is ~8 mag fainter at V (Section 2.3), implying that 
its contribution to combined light optical spectra is negligible. 
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Figure 1. Red optical spectrum of Ross 458AB (black line) com- 
pared to M0 (re d line), Ml (green line) and M2 (blue line) spectral 
templates from Bochanski ct al. (2007). All spectra are normal- 
ized at 7500 A, and templates are smoothed to the resolution of 
the observational data (A/AA 800) using a gaussian kernel. Re- 
gions of strong telluric absorption in the observational data are 
indicated. Inset box details the 6450-6800 A region highlighting 
the prominent Ha emission and absence of Li I absorption. 

A lower limit on the age of this system can be in- 
ferred from the optical spectrum of the primary. Fig- 
ure [1] displays A/AA ~ 800 optical spectroscopy for the 
combined AB components obtained by J. D. Kirkpatrick 
with the CTIO-1.5m RC Spectrograph on 22 May 199(£J 
The spectrum i s compared to M O , Ml and M2 spectral 
templates from iBochanski et al.l (|2007f l. the Ml show- 
ing the strongest agreement with the data. The promi- 
nent 6563 A Ha line is clearly seen in these data, but 
there is no indication of Li I absor ption at 6708 A, im- 
plying an age older than 30-50 M yr (Ma.ga zzu et al.lfl993t 
IChabrier et al.l 1996; Bildst en et al.lll997l ). Furthermore, 
there is no indication low surface gravity features, such 
as the weakened 7700 A K I doublet lines observed of 
coole r young M dwarfs up to the age of the Pleiades 
(e.g.. lMartin et al.lll996t Kirkpatrick et alll2008fl . These 
spectroscopic constraints place a lower limit on the age 
of this system of ^150 Myr. 

2.3. The Metallicity of the Ross 458 System 

The metallicity of Ross 458A can be determined 
from its position on the V — K/Mk color magni- 
tude diagram (CMD), a combination that effectively 
segregates mass and metalli c ity variations in field 
M dwarfs (iBonfils et all 120051: Uohnson fc Appsl 120091: 
ISchlaufman fc Laughlinl 120101) . For better accuracy, we 
disentangled the flux of Ross 458A from its closely- 
separated M7 com panion. Using AK = 4.41 from 
IBeuzit et all (|2004, a 0.1 mag uncertainty is assumed), 
we determined M Ka = 5.26±0.04 and M Kb = 9.67±0.11. 
We then used unresolved Hipparcos/2MASS photome- 
try (V - K = 4.180±0.016) and th e mean V - K/M K 
relation of ISchlaufman fc Laughlinl (j201Q[ ) to estimate 
AV" = 7.93±0.14 and hence (V - K ) A = 4.16±0.02 
The Mk /metallicity calibration of Uohnson fc Appsl 
(|2009t ) then yields [Fe/H] = +0.31±0.05, while the 

11 These data ar e available at DwarfArchives, 
|http: / /dwarf archives ■ org| 
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V — if /metallicity calibration of lSchlaufman fe Lau ghlin 
pOlCh yields [Fe/H] = +0.20±0.05. Hence, both rela- 
tions indicate a supersolar metallicity for the Ross 458 
system. However, these quantitative estimates should 
be regarded with care given the young age of the system. 
The primary may be system atically brighter /redder than 
the control samples used by I Johnson fe Appsl (|2009() and 
ISchlaufman fc Laughlin! (12010D . resulting in a systemat- 
ically higher metallicity. 

In summary, the Ross 458 system has a probable age 
of 150-800 Myr based on its rotation, strong activity 
indicators and spectroscopic features. It also appears to 
be metal-rich, with [Fe/H] = +0.2 to +0.3, albeit with 
possible age biases. We will compare these values to 
results from spectral modeling of the T dwarf companion 
in Section 4. 

3. NEAR-INFRARED SPECTROSCOPY 

3.1. The FIRE spectrograph 

FIRE is a cross-dispersed, near-infrared spectrometer 
recently installed at the auxiliary Nasmyth focus of the 
Magellan 6.5-meter Baade Telescope. It is a single-object 
spectrograph, designed primarily to deliver moderate- 
resolution (A/AA w 6000 for a O'.'6-wide slit) spectra 
over the 0.85-2.5 fim range in a single exposure. The 21 
orders dispersed by a combination of prisms and reflec- 
tion grating are imaged onto a 2048 x 2048, HAWAII-2RG 
(H2RG) focal plane array. A mirror can also be rotated 
in place of the reflection grating to obtain a single-order, 
high-throughput spectrum at lower resolution (A/AA = 
250-350; see below). A second H2RG camera images the 
entrance slit through an MKCFl J-band filter for source 
acquisition and manual guidi ng. Further details o n this 
instrument can be found in iSimcoe et al.l (|2008l I2010H 
and on the instrument webpageP^I 

3.2. Observations and Data Reduction 

Data presented here were obtained during FIRE's com- 
missioning run during 31 March through 7 April 2010 
(UT). Observations were obtained for three T dwarfs: 
Ross 458C, the T7.5 ULAS J141623.94+134836.3 
(aka SDSS J141624.08+134826 7B and here- 
after SPSS J1416 +1348B: IBurningham etH] 
[20101: ISchold I2010bt IBurgasser et al.l l2010bD and 

the T9 ULAS J 133553.45+113005. 2 [hereafter 

ULAS J1335+1130; IBurningham et al.l |200§ 1. A 
summary of the observations is given in Table [TJ 
Conditions ranged from clear (3 and 7 April) to patchy 
cirrus (6 April), with seeing of 0'.'6-0'.'8 at J-band. For 
all sources, we employed FIRE's prism-dispersed mode 
with the 0"6 wide slit aligned to the parallactic angle 
(rotator angle 89?5). The resolution for this setting 
ranges from A/AA w 350 in the 1.0-1.4 ^m Y- and 
J-band region to A/AA s» 250 in the 1.4-2.5 /im H- 
and -RT-band region, as determined from gaussian fits 
to emission lines in the arc lamp spectra and telluric 
emission lines in the science images. The spectrograph 
detector was read out using the 4-amplificr mode at 
"high gain" (1.2 counts/e - ). Each source was observed 

12 Maun a Kea Observatory fi lter s ystem; see Tokuna ga et all 
(2002) and Simons & Tokunaga (2002). 

13 http : //www . i irespectrograph. org 



in four 150 s exposures using an ABBA dither pattern, 
nodding 9" along the 30" slit, followed by observation 
of a nearby AO V star at similar airmass for telluric 
correction and flux calibration. Exposures of arc lamps 
(HeNeAr on 3 April; NeAr on 6 and 7 April) and quartz 
lamps, reflected off of the Baade secondary screen, were 
obtained with the target and AO V stellar observations 
for wavelength and pixel response calibration. 

Data were reduced using a combination of standard 
IRAF0 and custom IDlF°l routines. A pixel response 
frame was first constructed by combining two quartz 
lamp images taken with the lamp at high (2.5 V) and 
low (1.5 V) voltage settings. This was necessary to com- 
pensate for the substantial thermal emission from the 
lamp screen which tended to saturate the detector at the 
longest wavelengths. The two images were stitched to- 
gether at column 1300 (roughly 1.4 /im), then normalized 
by a median combination of all columns along the disper- 
sion region, which eliminated the intensity step function 
at the stitch and isolated pixel-to-pixel response varia- 
tions. The science images were divided by this pixel re- 
sponse frame, and bad pixels interpolated by averaging 
over nearest neighbors. The images were then pairwise- 
subtracted to remove sky background. Spectra were ex- 
tracted using the IRAF apall task. Calibrator stars were 
extracted first, then used as traces for the T dwarf and 
arclamp spectra. A wavelength solution was generated 
from a fifth-order spline fit to 60 lines in the arclamp 
spectra using the IRAF task identify, producing a mean 
dispersion of 8.5 A pixel -1 and scatter of 9 A («1 pixel). 
Note that the actual dispersion varies considerably along 
the chip, from 4 A pixel -1 at 1.0 fim to 20 A pixel -1 at 
2.4 /im. The individual spectra we re combined using th e 
xcombspec routine from SpeXtool (iCushing et al.l [20041) . 
using a robust weighted mean after normalizing in the 
1.2-1.4 /im region. The SpeXtool task xtellcor-general 
was used to compute telluric correction and flux calibra- 
tion from the AO V star, following the prescription of 
iVacca et al.l (|2003l ) and assuming a line profile width of 
18 A. 

To verify the flux calibration of our data, we com- 
puted MKO spectrophotometric colors from all three 
source s from the FIRE spectra following Cushin g"et al.l 
(2005), computing 100 iterations with each spectrum var- 
ied by its noise spectrum to determine uncertainties. 
For ULAS J1335+1130 and SDSS J1416+1348B, these 
colors match UKIDSS photometry to within the uncer- 
tainties, indicating robust relative flux calibration. For 
Ross 458C, the MKO colors are interme diate between 
the U KIDSS and ground-based GROND dGreiner et al.l 
120081) photometry reported in iGoldman et al.l (|2010l) , 
which themselves differ significantly (UKIDSS J — K 
= 0.51+0.06; GROND J - K = -0.21+0.06). This 
discrepancy may be due to differences in the filter 
basspands between UKIDSS, MKO and GROND fil- 
ter passbands, or may reflect intrinsic variability in the 
source. In Figure [2 we compare our reduced FIRE 

14 Image Reduction and Analysis Facility (IRAF; ITodvll 1989 1 is 
distributed by the National Optical Astronomy Observatory, which 
is operated by the Association of Universities for Research in As- 
tronomy, Inc., under cooperative agreement with the National Sci- 
ence Foundation. 

15 Interactive Data Language. 
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Figure 2. Spectral obse rvations of ULAS J1335+1130 based on Magcllan/FIRE (this study) and Gemini-North/NIRI plus VLT/ISAAC 
HBurnin gham ct al. 2008, red line). Both spectra are normalized at the 1.27 (im flux peak. The difference between these spectra (bottom 
panel) is less than 10% over most of the near-infrared band. 



Table 1 

Magellan/FIRE Observing Log. 



Source 


SpT 


UKIDSS J 


Date (UT) 


Airmass 


j- a 


AO V Calibrator 


Ross 458C 

ULAS J133553.45+113005.2 
ULAS J141623.94+134836.3 b 


T8 
T9 

T7.5p 


16.69±0.01 
17.90±0.01 
17.35±0.02 


3 Apr 2010 
7 Apr 2010 
6 Apr 2010 


1.34- 
1.38- 
1.40- 


-1.36 
-1.42 
-1.44 


642 
642 
642 


HD 114381 
HD 114381 
HD 124773 



a Total integration time divided into four ABBA dithered pointings. 

10.61 s frame readout. 

b aka SDSS J141624.08+134826.7B. 



spectrum of ULAS J13 35+1130 to that obtained by 
iBurningham et all (|2008l ) using the Gemi ni-North Near 
Infra red Imager and Spectrometer (NIRI; iHodapp et al.l 
|2003[ ) and Paranal Very Large Telesco pe Infrared Spec- 
trom eter and Array Camera (ISAAC; iMoorwood et al.l 
1998). The Burningham et al. spectrum was stitched to- 
gether from four orders individually flux-calibrated using 
UKIDSS YJHK photometry. Our spectrum is in excel- 
lent agreement with the NIRI and ISAAC data, with de- 
viations of less than 10% outside the telluric absorption 
regions. Importantly, the FIRE data reproduce the rel- 
ative YJHK flux peaks, which are essential for spectral 
model fits (Section 4). 



Total integration includes 



3.3. Results 



Figure [3] displays the reduced spectrum of Ross 458C, 
scaled to its absolute flux (F u units at 10 pc) based on 
UKIDSS J-band photometry and the Hipparcos parallax 
of the Ross 458 system. Signal-to-noise (S/N) in the 
spectrum ranges from 50 to 80 in each of the YJHK flux 
peaks. Ross 458C exhibits the classic spectral signatures 
of a T dwarf, with strong H2O and CH4 absorption bands 
at 1.15, 1.4, 1.6, 1.8 and 2.2 /mi; a steep red optical slope 
shaped by the pressure-broadened 0.77 /mi K I wings; 
and a relatively blue spectral energy distribution from 1.2 
to 2.1 /zm, du e in part to collision- induced H 2 absor ption 
(|Linskvl Il969t iSaumon etaLl H99l iBorvsowl [20021) . As 
anticipated from its extreme CHVfilter colors, the 1.6 /im 
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Figure 3. FIRE prism spectrum of Ross 458C in absolute flux units (F v ). The uncertainty spectrum is indicated by the grey line. 
Prominent H2O and CH4 absorption features are labeled, as well as regions influenced by the pressure-broadened K I doublet wing 
$ 1 A™ 1 ) an d collision-induced H2 opacity (A > 1.75 /im). The inset box shows continuum-normalized flux between 1.23 and 1.26 fim, 
highlighting the weak 1.24/1.25 /^m K I doublet present in the spectrum of this source (red lines show gaussian fits). 



CH4 band of Ross 458C is quite deep, while the J- and 
.ff-band flux peaks are narrow, all indicative of a very 
late-type, low-temperature T dwarf. Indeed, as shown in 
FigurelU these features in the spectrum of Ross 458C are 
intermediate between those of SDSS J1416+1348B and 
ULAS J1335+1130, indicating an intermediate spectral 
type. 

Our spectrum of Ross 458C also reveals the pres- 
ence of the 1.24/1.25 fim K I doublet, highlighted in 
the inset box in Figure [3] Th ese lines are known to 
fade in the latest-typ e T dwarfs ()Burgasser et al.i r2002a; 
IMcLean et al.l 120031) . likely a consequence o f KCl(g) 
forma tion at temperatures below 1000 K (jLoddersI 
jl999[). We measured pseudoequivalent widths pEW = 
2.7±0.2 A for each doublet line, comparable to measure- 
ments for other T7 -T8 dwarfs ([Burgas ser et al.l r2002at 
IMcLean et al.l [2003h . We do not see these lines in our 
spectra of ULAS J1335+1130 or SDSS J1416+1348B 
to limits of 0.5 and 0.3 A, respectively. For the lat- 
ter source, this may be a consequence of its subso- 
lar metallicity in analogy to t he metal poor T 6 dwarf 
2MASS J09373487 +2931409 (IBureasser et al.l [2502at 
IMcLean et "aLll20Q7t see Section 4.2) 

Another feature of note in the spectrum of Ross 458C 
is its fairly pronounced 2.1 /im i^-band peak. The rel- 
ative brightness of this peak is greater than those of 



both ULAS J1335+1130 and SDSS J1416+1348B, in ac- 
cord with its reddest J — K color. The i^-band peak 
is normally suppressed by H2 absorption in late-type T 
dwarfs. This absorption is weakened in brown dwarfs 
with low surface gravities (i.e., young, low-mass) and/or 
higher metallicities, as both effects result in a lower pho- 
tospheric pressure and weaker collision-induced absorp- 
tion. For Ross 458C, both of these properties likely con- 
tribute based on the characteristics of the system (Sec- 
tion 2). 

3.4. Spectral Classification 

To classify the T dwarfs, we used the near-infrared 
spectral indices and index-spec tr al type relati o ns de - 
fine d in IBureasser et all (I2006blh IWarren et all (120071) : 
and IBurningham et all (|2008f ) : these are listed in Ta- 
bled Our index measurements for ULAS J1335+1130 
and SDSS J1416+13 48B are in rough agr e ement 
with previous r esults (IBurningham etaU 120081 [2010; 
IBureasser et~aH l2010bl) . with minor variances at- 
tributable to differences in spectral resolution and S/N 
between data sets. We derive mean classifications of 
T8, T8.5 and T7.5p for Ross 458C, ULAS J1335+1130 
and SDSS J1416+1348B, respectively. Uncertainties of 
0.5 subtypes are assumed for the first two sources, while 
the last classification is deemed "peculiar" due to the 
larger scatter in indices (±1.2 subtypes). The last two 
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Figure 4. FIRE prism spectra of (top to bottom) Ross 458C, ULAS J1335+1130 and SDSS J1416+1348B. Data are normalized at 1.27 
and offset by constants (dotted lines). Note the variations in the Y- and i-f-band flux peaks. 



Table 2 

Spectral Indices for Observed T Dwarfs. 



Index 


Ross 458C 


ULAS J1335- 


rll30 


SDSS J1416+1348B 


Ref. 


Value 11 


SpT 


Value- 1 


SpT 


Value" 


SpT 


H2O-J 


0.051±0.002 


T8.0 


0.028±0.003 


T8.5 


0.044±0.003 


T8.0 


1 


CH4-J 


0.154±0.001 


T8.5 


0.109±0.002 


T9.0 


0.262±0.002 


T7.0 


1 


Wj 


0.272±0.002 


T9.0 


0.253±0.002 


T9.0 


0.360±0.002 


T7.0 


2.3 


H 2 0-H 


0.193±0.003 


T7.5 


0.192±0.007 


T7.5 


0.148±0.004 


T9.0 


1 


CH4-H 


0.112±0.003 


T8.0 


0.068±0.005 


T9.0 


0.192±0.003 


T7.0 


1 


NH3-H 


0.632±0.004 




0.537±0.007 




0.610±0.005 




1 


CH4-K 


0.078±0.005 


T7.5 


0.073±0.013 


T7.5 


0.191±0.032 


T5.5 


1 


K/J 


0.167±0.001 




0.130±0.002 




0.044±0.001 




1 


J - H b 


-0.14±0.01 




-0.34±0.02 




-0.25±0.01 






H — K b 


0.38±0.02 




0.06±0.05 




-1.02±0.09 






J - K b 


0.24±0.02 




-0.28±0.04 




-1.25±0.09 







References. - (1) Burgassc r et all J200 6b1; (2) IWarren et all {2007); (3) 
IBurningham et al.l l l200g) ; f4) IDelorme et al.l 1 12008) . 



a Values for each source were measured for 1000 realizations of the spectrum (100 
realizations for spectrophotometric colors), each with a normal distribution of ran- 
dom values scaled by the noise spectrum added to the original fluxes. The reported 
values are the means and standard deviations of the se measurements. 
b Spectrophotometric colors computed according to Cushing et al. (2005). 
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classificat ions formally agree with previously published 
results (|Burningham et all 120081 120101 IBurgasse r et al.l 
l2010bft . Our classification of Ross 458C splits t he dif- 
ference betw een the photometric estimates of iScholzl 
(|2010aft and iGoldman etafl (|2010t l. and is consistent 
with the relative ordering of the spectra in Figure 21 The 
absolute J magnitu de of Ross 458C is also comparable 
to other T8 dwarfs (|Vrba et al.ll200l . 

3.5. Physical Properties of Ross 458C 

With a spectral type and known distance, it is pos- 
sible to assess some of the basic physical properties of 
Ross 458C. A bolometric luminosity was determined 
using the absolute UKIDSS JHK- b&nd magnitude s 
and bolometric corrections (BC) from lLiu et al.l (12010). 
Adopting BCj = 2.62±0.14, BC H = 2.19±0.07 and 
BCk = 2.05±0.09 for a T8 dwarf, we derive estimates of 
log 10 L bol /L Q = -5.69±0.06, log 10 L bol /L Q = -5.65±0.04 
and log 10 Lboi/Lo = -5.55±0.05 based on J-, H- and 
K-hand fluxes, or a weighted mean of log 10 L bo i / X© = 
-5.62±0.03. This luminosity is similar to that of the 
T8 spectral standard 2MASS J04151954-0935066 (here- 
after 2MASS J0415 -0935; log 10 L bol /L Q = -5.67±0.02; 
iSaumon et al1l2007ft . The T e// of Ross 458C can be 
inferred from its luminosity by adopting a radius from 
the evolutionary models of SM08. For an age of 150- 
800 Myr (Section 2.2), these models predict a radius of 
0.12-0.13 R©, T e ff = 650±25 K and a mass in the range 
0.006-0.011 M©. The inferred temperature is cooler 
than that of 2MAS S J0415-0935 (T e// = 750±25 K; 
ISaumon et al.l 12007ft due to the inflated radius, while 
the inferred mass is below th e deuterium burning min- 
imum mass limit (0.014 M©; iChabrier fc Baraffd 120001 : 
ISpieeel et al1l20loh . 

4. SPECTRAL MODEL FITS 

4.1. Spectral Models and Method 

We compared the FIRE spectra of these three 
T dwarfs to a s uite of atmosphere models from 
IMarlev et aLI (l200l and SM08 follow i ng the prescrip- 
tions d etaile d in ICushing et al.l (|2008ft : IBurgasser et al.l 
(|2008bft and lBowler et al l (|2009ft . The one-dimensional 
models include all sig n ificant gas species consid- 
ered in iFreedman et all (|2008fl as well as C0 2 (cf. 
lYamamura et al.l l2010f ) . albeit with line lists for CH4 
and NH 3 that are known to be incomple te (e.g., 
ISaumon et al.l l2006h . Solar abundances from lLoddersI 
(2003) were assumed. Non-equilibrium chemistry of CO, 
CH4, H2O and NH3 (the relevant molecules influencing 
the near-infrared spectra of T dwarfs) by vertical trans- 
port was accou nted for using an eddy diffusion parame- 
teriza t ion. K„ (IFeeley fc Lodderslll99l IGriffith fc Yelld 
[l999t ISaumon et all 120061 : iHubenv fc Burrows! 12007). 
Equilibrium condensate formation was included in 
the chemistrjO, with grain size and vertical distri- 
butions parameterized by a sedimentation rate f sec i 

16 Error estimates quoted here take into account uncertainties in 
the spectral classification (±0.5 subtype), photometry, parallax of 
Ross 458 and the spectral type/BC relations (0.14, 0.07 and 0.08 
mag for J, H and K, respectively). Systematic errors — e.g., from 
gravity or metallicity effects — are not accounted for. 

17 The chemical equilib rium calculation accounts for all im- 
porta nt condensed species JLodders fc Fcglcv 2006; Visschc r et al.l 
12010 ). For the cloud opacity we included only silicates, iron, and 



(|Ackerman fc Marlevi |200H: IMarlev et all 12002ft . While 
contemporaneous studies of late-type T dwarfs have 
gener ally ignored the opa c ity fr o m condensate grains 
/e.g.. iLeggett et all 120071 120091 : ISaumon et all 120071 : 



IBurgasser et al.M2008bf ). we considered two sets of SM08 
models with ("cloudy") and without ("cloudless") con- 
densate opacity included. Both sets of models use the 
same gas chemistry and have the same constituent abun- 
dances at any given pressure and temperature. The 
cloudless models assume that condensates do not con- 
tribute to the opacity (but do contribute to the gas chem- 
istry) , perhaps because they have "rained out" of the at- 
mosphere or occupy only a small fraction of the visible 
surface area. The cloudy models explicitly include opac- 
ity from cloud grains. The cloudless models sampled at- 
mospheric parameters T e ff = 500-1000 K (50 K steps); 
log g = 4.0-5.5 cgs (0.5 dex steps); [M/H] = -0.3, and 
+0.3 dex relative to solar; and K zz = and 10 4 cm 2 s _1 . 
The cloudy models sampled the same parameter space 
with f sec i = 2 and solar metallicity alone. Hereafter, we 
refer to sets of parameters by the notation k = (T e tf, 
logs, [M/H], f sed , logK zz ), where f sed = NC for the 
cloudless models. The model spectra were computed as 
surface fluxes (F v flux units), smoothed to the mean res- 
olution of the FIRE prism data (A/ A A = 300) using a 
gaussian kernel, and interpolated onto the same wave- 
length scale as the data. We also assigned physical pa- 
rameters of mass, age and radius to each spectral model, 
interpolating the appropriate evolutionary tracks from 
SM08. 

The observational data (also in F v flux units) were 
scaled to their apparent UKIDSS J-band magnitudes. 
Simultaneous fits were made to the 0.85-1.35 /im, 1.42- 
1.8 /im and 1.95-2.35 /im regions to avoid strong tel- 
luric absorption. For each data-model compar ison, the 
goodness-of-fit statistic Gk (|Cushing et al .112008ft was cal- 
culated following the same weighting scheme employed 
in that study, with each pixel weighted by its breadth 
in wavelength space (wj cx AXj). Model surface fluxes 
were scaled by the factor Gk = (R/d) 2 which minimizes 
Gk (Equation 2 in ICushing et al.ll2008ft . where R is the 
radius of the brown dwarf and d its distance from the 
Earth. As in IBurgasser et all (|2008bft . we generated dis- 
tributions of the fit parameters by weighting each model's 
contribution according to the F-distribution probability 
distribution function (PDF): 



Wi oc 1 - F(G ki /min(G k ) | v eff ,v eff ). 



(1) 



Here, Gki/min(Gk) is the ratio of goodness-of-fit (effec- 
tively a ratio of chi-square residuals) between the best-fit 
model and the i th model, and F(x\u e ff,u e ff) is the F- 
distribution PDF for effective degrees of freedom 



max({ra}) 



N 

E 

4=1 



(2) 



The sum is oyer all pixels included in the fit (see 
IBurgasser et alll2010aft . The F-distribution provides a 

AI2O3 (alumina or corundum). Other condensates are expected 
to form above the silicate cloud at low effective temperature, but 
of t hese only Na2 S is likely to have any significant column opac- 
ity (Marlcy 2000). In a future work we will consider the possible 
contribution of this species. 
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more robust assessment of the equivalence of different 
model fits than t he exp onenti al weighting em p loyed i n 
iBurgasser et~all (|2008b[ ) and iBurgasser etail (|2010b[ ). 
Parameter means and one-sided standard deviations were 
separately computed as 



(k) 



EiWi 



and 



E t± ^(k,-(k)) 2 



(3) 



(4) 



The first sum was performed over all models and the 
second restricted to those models whose parameter values 
are above (+) or below (— ) the weighted mean (k). 

We examined several variations in our fitting method 
to assess internal systematic effects. Modest changes to 
the spectral regions that contributed to computing Gk 
made little difference to the overall best fit. For ex- 
ample, including the telluric regions in the fits, or ex- 
panding the lower and upper ends of the spectral range 
used, resulted in changes of less than 10 K and 0.1 dex 
in the mean T e ff or \ogg values for all three sources. 
Only by excluding whole spectral regions — such as the 
individual YJH A'-band flux peaks — did inferred T e ff 
and log g values change significantly, by up to 100 K and 
0.5 dex (similar to band-to-band variations reported in 
iCushing et al.l 120081 ). Variations in the pixel weighting 
scheme also had negligible effect. Weighting all pixels 
equally gave identical parameters (to within 10 K and 
0.1 dex in T e ff and log g) as weighting by pixel bandpass. 
Finally, using the exponential parameter weighting em- 
ployed in IBurgasser et al.l (|2010bl ). Wi oc e - n - 5G k i; also 
produced negligible changes in the inferred parameters 
and their uncertainties. We therefore assume systematic 
uncertainties of 20 K in T e ff and 0.2 dex in logg, al- 
though we cannot rule out larger systematic biases based 
on the accuracy of the models. 

4.2. Results 

Table [3] summarizes the inferred parameters for each 
source, with results for the cloudless and cloudy models 
reported separately. Figures [SHU show the best-fit mod- 
els to the individual spectra, while Figures I5HTU1 show 
the parameter distributions. It is first worth noting the 
quality of the fits. For both cloudless and cloudy models, 
deviations are seen in regions of strong CH4 absorption, 
particularly in the 1.6-1.7 /im region for Ross 458C and 
ULAS J1335+1130 and the wings of the 1.27 ^m J-band 
peak for SDSS J1416+1348B. As noted previously, these 
discrepan cies are likely the r esult of incomplete CH4 
opacities (Saumon et al. 2006). However, there are also 
differences between the cloudless and cloudy model fits, 
notably in the relative YJH peak fluxes of Ross 458C 
and ULAS J1335+1130. In the cloudless model fits of 
these sources, the -ff-band peak appears to be suppressed 
relative to J, while the V-band peak appears suppressed 
in the cloudy model fits. These differences are less pro- 
nounced in fits to SDSS J1416+1348B. The cloudless 
models also exhibit excessively strong K I lines at 1.24- 
1.25 fim, while these lines in the cloudy models are more 
consistent with the data. For both models, the 2.1 /im 
A-band peak is reasonably well-reproduced for all three 



sources (see also ICushing et all 120081 ). Nevertheless, it 
appears that the cloudy models provide better matches 
to the data than the cloudless ones. 

Quantitative differences between the cloudless and 
cloudy model fits are also seen. For Ross 458C (Fig- 
ures [5] and [5]), the cloudy model provides a statisti- 
cally superior fit to the data, significant at the 99.99% 
confidence level based on the F-test PDF. The cloudy 
model also predicts a significantly cooler temperature 
than the cloudless case, T e ff — 6351^ K versus T e ff 
= 760^4° K. The former value is consistent with the 
luminosity-based T e ff determined in Section 3.5, while 
the latter is consistent with a "normal" T8 field dwarf 
like 2MASS J0415- 0935 (e.g., IGolimowski et~aTl l200i 
Steph ens et al.ll2~009f) . Inferred surface gravities are com- 
parable for the two models, favoring log g « 4.0 cgs, while 
the cloudless models indicate a metal-rich atmosphere 
([M/H] pa +0.3). These parameters are consistent with 
the youth and metallicity of the Ross 458 system, as the 
mean cloudless and cloudy model ages are 0.3toJ Gyr 
and 0.14^01 Gyr, respectively. Both models also predict 



+0.006 
-0.004 



M and 



a very low mass for Ross 458C, M = 0.010 
M <0.006 M Q , again at or below the deuterium burning 
mass limit. There is a modest difference between whether 
vertical mixing is slightly favored (cloudless models) or 
not (cloudy models), but it is not statistically signifi- 
cant. The cloudless and cloudy models predict distinct 
distances of 15.6±2.5 pc and 9.4±0.7 pc, respectively, 
based on the optimal model scaling. These differ by 2.6cx, 
and are 1.8er higher and 2.4cr lower than the parallactic 
distance of the Ross 458 system, although the latter is 
off by only 11%. 

For ULAS J1335+1130 (Figures M and Mj, the best-fit 
cloudy model also has a smaller Gk value than the best-fit 
cloudless one, although the statistical significance is re- 
duced (84% confidence level). Again, the cloudy models 
predict a significantly lower effective temperature than 
the cloudless models (T e ff = 595l 4 g K versus T e ff = 



700 



+30 
-50 



K). The former is consistent with previously re- 



port e^_j^alues_by_ 



and lLeeeett et al 



Burningham et al.l (pOOl 550-600 K) 



(2009, 



500-550 K). The cloudy mod- 
els also indicate a lower surface gravity (logg « 4.0 cgs 
versus logg = 4.7±0.4 cgs), while differences in vertical 
mixing are statistically unconstrained. Like Ross 458C, 
the cloudless models converge on a supersolar metallicity 
for ULAS J1335+1130, [M/H] = +0.23± 1 ;^. The differ- 
ing parameters between the cloudless and cloudy models 
translate into marginally overlapping ages (2.0+^g Gyr 
versus 0.18±j^ Gyr) and masses (0.010-0.032 M ver- 
sus 0.005-0.013 M ), and distinct distances (20.5±3.6 pc 
versus 14.4±1.1 pc, a 1.6(7 discrepancy). 

Finally, for SDSS J1416+1348B (Figures and [T0j), 
we again find that the best-fit cloudy model is marginally 
better than the best-fit cloudless one, albeit at even lower 
statistical significance (confidence level of 77%) reflecting 
the similarity in the near-infrared SEDs of these models. 
Nevertheless, the inferred temperatures are marginally 
distinct, with the cloudy models again indicating a cooler 
temperature (T e// = 595±|g K versus T e// = 685±|| K). 
Unlike Ross 458C and ULAS J1335+1130, the inferred 
surface gravity of this source is considerably higher (log g 
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Figure 5. Best- fit cloudless (top) and cloudy (bottom) model fits (red lines) to the observed spectrum for Ross 458C (black lines). 
Observed data are scaled to the apparent UKIDSS J-band magnitude of the source, while the model is scaled to minimize Gk (indicated, 
along with d/R flux scaling factor). The noise spectrum for Ross 458C is shown by the grey lines, and fit regions indicated by the shaded 
boxes at the top of each panel. 
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Figure 6. Same as Figure[5]for ULAS J1335+1130. 
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Figure 7. Same as Figure[5]for SDSS J1416+1348B. 



Clouds in Ross 458C 



13 



A 0.8 - 



S 0.8 - 



500 600 700 800 900 1000 
T.„ (K) 



4.5 5.0 5.5 
log g (cgs) 



3 10 a 

K zz {cm 2 s' 1 ) 



-0.3 +0.3 
[M/H] (dex) 



1.2r 

c 1.0- 
o 

1 0.8- 

N 

| 0.4- 
o 

2 0.2- 
0.0 L 



500 600 700 800 900 1000 
T e ff (K) 



4.5 5.0 5.5 
log g (cgs) 



10" 

K zz {cm 2 s' 1 ) 



Figure 8. Parameter distributions for cloudless (top) and cloudy (bottom) model fits of Ross 458C, based on the parameter weighting 
scheme described in the text. From left to right, distributions in T e jy, log g and K zz are shown. For the cloudless models, the distribution 
in [M/H] values is also shown. 
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Figure 9. Same as Figure|8]for ULAS J1335+1130. 
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Figure 10. Same as Figure|8]for SDSS J1416+1348B. 



= 5.2±0.3 cgs cloudless; log g ~ 5.5 cgs cloudy), while the 
cloudless model indicates a subsolar metallicity ([M/H] 
= - O.Uto H). As di s cussed in iBurningham et al.l (|2010H 
and lBurgasser et al.l (|2010bD . high surface gravity and/or 
subsolar metallicity are required to produce both the sup- 
pressed -fT-band peak and broadened V-band peak ob- 
served in the spectrum of this source. Vertical mixing 
is again poorly constrained. The inferred atmospheric 
parameters of SDSS J1416+1348B indicate a brown 
dwarf that is older (4-12 Gyr) and more massive (0.03- 
0.05 M Q ) than the other two sources, consistent with 
prior analyses. Surprisingly, the cloudless and cloudy 
models converge on a common distance of 11±3 pc, larger 
than but formally c onsistent with t he preliminary paral- 
lax measurement of lSchofzl (|2010bL 7.9±1.7 pc). 

5. DISCUSSION 

5.1. Clouds in Young, Cold and Metal-Rich T Dwarfs 

The above analysis illustrates how condensate opac- 
ity, when included in spectral models, can have a pro- 
found effect on the near-infrared spectral energy dis- 
tributions and inferred T e //s of the coldest T dwarfs. 
More importantly, the inclusion of condensate opacity 
provides a statistically significant improvement in fits to 
the near-infrared spectrum of Ross 458C, and marginally 
so for ULAS J1335+1130 and SDSS J1416+1348B. This 
is the not the first indication of condensates being present 
in T dwarf atmospheres. Early-ty pe T dwarfs exhibit 
evidence of waning cloud opacity (iMarlev et all 120021 : 
iCushing et al.ll2008t ; ILe"ggett et all2008lh and SM08 have 
shown that T dwarf models produce excessively blue 
near-infrared colors without the inclusion of condensate 
extinction. Nevertheless, several theoretical studies have 
concluded that condensate clouds play little or no role in 
shaping the near-infrared spectra of brown dwarfs cooler 
than wlOQQ K, as they reside below the visible p hoto- 
sphere fe.g.. lAllard et al.ll200lt iBurrows et al.ll2006D . In- 



stead, we find that clouds are essential for reproducing 
the near-infrared SED of the T e// w 650 K T dwarf 
Ross 458C, and to a lesser degree the T e ff ~ 600 K T 
dwarfs ULAS J1335+1130 and SDSS J1416+1348B. 

The importance of clouds in these two sources is 
likely related to their youth and associated low sur- 
face gravities. Among the L dwarfs, a correlation 
between low surface gravity and thi ck clouds has 
been n oted by ICushing et al.l (|2008l ) and iStephens et al.l 
(2009 ) based on broad-band s pect ral model fits: and by 
iMetchev fc Hillenbrand! (|2006[ ) and lDupuv et alj ()2009bl ) 
based on the delayed disruption of clouds at the L 
dwarf/T dwarf transition. There is also evidence that 
this correlation continues into the T dwarf regime, based 
on cloudy model fits to broad-band spectral data for 
the low-gravity T5. 5 d warf SDSSp J111010.01+011613.1 
(jStephens et al.l |2009[ ). Finally, the 9-11 ^m silicate 
grain feature has also been shown to be more pro- 
nou nced in L dwarfs with unusually red J — K col- 
ors ([Burgasser et al.ll2008at lLooper et al.ll2008| ). sources 
whose kinematics reflect youth and low s urface gravities 
(|Fahertv et al.ll2009t [Schmidt et al.ll2010D . 

Atmospheric models show that two complementary ef- 
fects can drive the persistence of clouds in low gravity, 
low-temperature brown dwarf atmospheres. First, the 
column of atmosphere above a given pressure level is 
deeper in a lower gravity dwarf, and as a consequence 
the atmosphere is everywhere warmer at a fixed T e ff- 
This implies that silicate clouds lie higher in these at- 
mospheres, providing a larger column depth of extinc- 
tion in regions of mini mum molecular gas opacity. Sec- 
ond, cloud m o dels f rom Ac kerman fc Marleyl ([200 1| ) and 
iHelling et al. (2008!) both predict that mean cloud par- 
ticle sizes are larger in lower gravity atmospheres. In a 
T e f f = 650 K SM08 model, for example, the mean par- 
ticle radius is roughly ten times larger for log g = 4.5 
as compared to log g = 5.5. The sub- micron particles 
found in higher gravity models contribute little opac- 
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Table 3 

Model-fit Parameters for Observed T Dwarfs. 



Parameter Best Fit Mean Best Fit Mean Prior Value Ref. 

(Cloudless) (Cloudless) (Cloudy) (Cloudy) 



Ross 458C TT8J 





71.8 




51.3 








d/R (pc/R Jtl p) 


11.5 


13.4±1.9 


7.9 


7.5±0.6 






Te// (K) 
logg (cm s~ 2 ) 
[M/H] (dex) 
K zz (cm 2 s" 1 ) 


700 
4.0 
+0.3 
10 4 


760+™ 
4 2+°' 3 
fs+0.3 


650 
1.0 

a 



635+J* 
«4.0 
a 


+0.2. ..+0.3 


1 


Mass (M ) 


0.006 


01 0+ 006 
u.uiu_ 004 

n + 0.5 
U.d_ 2 


0.006 


<0.006 


0.006-0.011 b 


1 


Age (Gyr) 


0.08 


0.13 


14+ 003 


0.4-0.8 
0.4-2 


1 
2 


R (Rjup) 

d (pc) 


1.22 
14.1 


1 16+ ' 07 
J - 1D -0.09 

15.6±2.5 


1.25 
9.8 


ssl.24 
9.4±0.7 


11.7±0.2 


3 


ULAS J1335+1130 (T9) 




22.7 




21.0 








d/R (pc/R Jup ) 


20.9 


20.2±2.4 


12.0 


11.8±0.8 






Te// (K) 


700 


rao±l° 


600 




500-600 


4,5 


log g (cm s 2 ) 


5.0 


4.7±0.4 


1.0 


«4.0 


4.0-4.5 


1,5 


[M/H] (dex) 
K zz (cm 2 s" 1 ) 


+0.3 
10 4 


+0.23lg;g 


a 



Q a 


0...+0.3 


5 


Mass (M©) 


0.031 


0.021±0.011 


0.006 


0.005-0.013 


0.014-0.030 


4 










0.005-0.019 


5 


Age (Gyr) 


3.3 


2.0t 2 ;j? 


0.17 


n io+0.23 


0.6-5.3 
0.01-2 


1 

5 


R (Rjup) 


0.90 


1 01+ ' 13 
i - ul -0.12 


1.23 


fsl.23 


1.07-1.17 


5 


d (pc) 


18.7 


20.5±3.6 


14.7 


14.4±1.1 


8-12 


4,5 






SDSS J14164 


-1348B (T7.5p) 








22.3 




21.0 








d/ii (pc/R Jup ) 


15.0 


13.2±3.6 


15.1 


13.2±4.0 


12.8±3.0 


6 


Te// (K) 


700 


685l« 


600 




650±60 
500 


6 
7 


logg (cm s~ 2 ) 


5.5 


5.2±0.3 


5.5 


w5.5 


5.2±0.4 


6,7 


[M/H] (dex) 
K zz (cm 2 s" 1 ) 


0.0 

10 4 


17+ - 17 
-"■'-1-0.13 


a 





-0.3 
10 4 


6,7 
6,7 


Mass (Mq) 


0.052 


o n4i+°- 011 


0.036 


noc+O.OOl 

u.uoo_ 003 


0.021-0.045 
0.029-0.038 


6 
7 


Age (Gyr) 


12 


8±4 


12 


6-12 


2-10 


6.7 


R (Rjup) 


0.78 


0.84±0.06 


0.86 


«0.86 




6 


d (pc) 


11.8 


11.1±3.2 


13.0 


11.4±3.4 


7.9±1.7 


8 



Refere nces. — (1) This p aper; (2)lGoldman et all |2010l ) : (3)lvan Leeuwenl ll2007[) : 
(4) iBurningham et all jl2008j) : (5) iLeggett et all l|2009l ): (6) iBurgasser et alJ (|2010bl ): 
(7) Bumiiigham et~aT1 IpOUl ); (8 HSchofzl (|2010bD . 
a Solar metallicity is assumed for all cloudy models. 

b Mass estimate based on the luminosity and estimated age of Ross 458C, 
cloudy evolutionary models of SM08. 
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ity in the near-infrared because of Mie-scattering effects, 
rendering the deep silicate cloud less visible. Clouds in 
lower gravity dwarfs are therefore more visible and scat- 
ter more efficiently than those in their older and more 
massive counterparts. This interpretation explains why 
cloudless models have so far provide d adequate fits to the 
spectra of log a > 5 T dwarfs (e.g.. iGeballe et al.M2001t 
iSaumon et al.l 2007; B urgasser et alJ l2008bf ). including 
SDSS J1416+1348B. 

Metallicity is another factor to consider. Cloud for- 
mation appears to be inhibited in metal-poor L sub- 
dwarfs, as exemplified by the retentio n of Ca, Ti and 
TiO gas to cooler TV. f f s in these ob jects (jBurgasser et al.1 
\m)3. [20071 IGizis fc Harvinl 120061 Similarly, unusu- 
ally blue L dwarfs — including the brighter compan- 
ion of SDSS J1416+1348B — appear to have slightly 
metal-poor atmospheres with unusually thin clouds 



(jBurgasser et al.1 l2008al iCushing et al.l I2010D or frac- 
tured cloud layers (jMarlev et al.l IsubmittedT ) . It is not 
unreasonable to speculate tha t supersolar metal licity 
dwarfs may have thicker clouds (jLooper et al.ll2008h , yet 
theoretical predictions are mixed. SM08 have shown 
that warmer cloudy models (T e ff > 1300 K) appear 
"cloudier" (redder) for supersolar metal licities. How- 
ever, c hemical equilibrium calculations by iVisscher et al.1 
(2010) find that MgSiOa (enstatite) condenses at warmer 
temperatures in higher metalllicity atmospheres, leading 
to a deeper cloud that may be buried by gas opacity 
in cooler brown dwarf s. Atmospheric model calculations 
bv lWitte et al.l (|2009t ) find that grain size and cloud op- 
tical depth do not change appreciably between [M/H] = 
and +0.5 for low-temperature brown dwarfs, although 
the clouds themselves expand into lower pressure regions 
of the atmosphere and may have a greater effect at wave- 
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lengths corresponding to minimum gas opacity. Metal- 
licity effects on cloud formation and cloud properties in 
cool T dwarf atmospheres clearly remains an area of fur- 
ther theoretical exploration. 

In summary, the unexpected presence of condensate 
absorption in the spectrum of Ross 458C (and possi- 
bly ULAS J1335+1130) can be attributed to its low 
surface gravity, and possibly its supersolar metallic- 
ity. This hypothesis follows similar trends observed in 
the L dwarfs. It also explains why clouds had not 
previously been needed in model fits of high surface 
gravity T dwarfs. The importance of clouds in shap- 
ing low gravity T dwarf spectra should nevertheless 
be a consideration in the characterization of planetary- 
mass brown dwarf candidates in local young clusters 
(e.g.. iZapatero Osorio et alJ l2002t iCasewell et all 120071 : 
iBurgess et al.1 120091 iMarsh et al.1 120101 ) and low -mass 
cxopl anets recently imaged around nearby young stars 
(e.g-. lKalas et alj|2008t iMarois et a l.ll2008t IJa nson et al.l 
MMk lLagrange et al.ll2oToli ~ 

5.2. Clouds and the Convergence of Near-infrared and 
Mid-infrared T Dwarf Temperature Scales 

Spitzer has enabled detailed studies of cold T dwarfs in 
the mid-infrared, wh ere the majority of the SED emerges 
(|Leggett et al.ll2010l ). Model fits to mid- infrared spectral 
data have proven particularly robust, with agreement 
in continuu m fluxes and prominent H2O, CH4 and NH3 
bands (e.g.. iSaumon et al.l 120071 : IBurgasser et al.ll2008tt 
ILeggett et al.1 120091 ) . However, there have been indica- 
tions that temperatures inferred from near- and mid- 
infrared analyses can d iffer significantly . A case in point 
is SDSS J1416+1348B. IBurgasser et all (|2010bD inferred 
T e ff = 650±60 K for this source from near- infrared 
cloudless spectral mo del fits, similar to the clou dless re- 
sults presented here. iBurningham et all (|2010l ) inferred 
T e ff ~ 500 K based on the extreme H — [4.51 co lor of 
this source ([Warren et alj|2007t ILeggett et al.ll2010f ). We 
find that inclusion of cloud opacity in the near-infrared 
modeling of this source brings inferred T e ffS closer to 
mid-infrared results, despite having minimal effect on the 
shape of its spectrum. Even better convergence is seen 
for ULAS J1335+1130, as our near-infrared cloudy model 
temperature is f ully consistent with estim ates based on 
H — [4.5] color (|Burningham et all 120081) and near- t o 
mid-infrared spectral model fits ([Leggett et all [2009). 
For Ross 458C, the cloudy model temperature is consis- 
tent with that inferred from evolutionary models based 
on the age and luminosity of the source (Section 3.5). 

Prior mid-infrared analyses have been based on cloud- 
less model fits, however, so does the inclusion of clouds 
affect these results? To explore this question, Fig- 
ure QT] compares the observed mid-infrared spectrum of 
ULAS J1335+1130 o btained with the Spitzer Infrared 
Spect rograph (IRS; iHouck et all 120041: ILeggett et all 
2009) to three SM08 models: two sampling the best-fit 
parameters from our near-infrared cloudless and cloudy 
model fits, and one sampling the c loudless model fit 
parameters from Leggett et al. (2009)0 The first two 
models are scaled according to the near-infrared nor- 

18 ILeggett etahl (20031) used K zz = 10 6 cm 2 s _1 in their best- 
fit models whereas we use K zz = 10 4 cm 2 s — 1 in Figure [TT] The 
difference between these models is negligible. 
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Figure 11. Spitze r /IRS spectrum of ULAS J1335+1130 from 
ILeggett et al.l (20091 . black line, n oise spectrum in grey), scaled 
to its apparent 8.0 (im magnitude (Burningham ct al. 2008), com- 
pared to best-fit cloudless (blue line) and cloudy (solid red line) 
models from our near-infrared analysis, and the be st-fit cloudless 
model from the near- and mid-infrared analysis of Legg ett et al.l 
(2009, green line). Model surface fluxes have been smoothed to 
the resolution of the IRS data (A/ A A PS 100) using a Gaussian 
kernel, and scaled according to their respective best-fit distances 
and radii. The best-fit cloudy model to the near-infrared data with 
K zz = 10 4 cm 2 s _1 is also shown (dashed red line). Observational 
uncertainty in flux at 10 (im in indicated by the error bar. 

malizations from our analysis, while the Leggett et al. 
model is scaled to the best-fit radius and distance de- 
rived in that study. The 600 K cloudy model comes much 
closer to reproducing the pseudo-continuum flux of the 
observed data than the 700 K cloudless model, but is still 
underluminous and predicts excessively strong NH3 ab- 
sorption in the 9-13 fim region. This is not a cloud effect, 
however, but the absence of vertical mixing in the best- 
fit cloudy model. Vertical mixing reduces NH3 and CO 
photospheric abundances through non-equilibrium chem- 
istry, but has little i mpact on the near-inf r ared spectrum 
for T eff < 1000 K ([Griffith fc Yelld[T999HSaumon et al.1 
2006). This is confirmed in our spectral modeling, as the 
presence or absence of vertical mixing is statistically am- 
biguous (Figures l8l- [T0]) . It has a profound effect on mid- 
infrared spectra, however, as the T e ff = 600 K, logg = 
4.0, K zz = 10 4 cm 2 s" 1 cloudy model is a good match to 
data for ULAS J1335+1130 (i.e., consistent within the 
observational uncertainties). As such, the inclusion of 
cloud opacity does not appear to modify mid-infrared 
analyses for the coldest T dwarfs. 

Convergence between near-infrared and mid-infrared 
spectral model fit parameters and temperature scales for 
cold T dwarfs is important for both near- and long- 
term study of these objects. With Spitzer now de- 
pleted of cryogen, our ability to obtain mid-infrared 
spectroscopy for faint brown dwarfs is restricted pend- 
ing deployment of the James Webb Sp ace Telescope 
(JWS T) Mid-Infrared Instrument (MIRI; iWright et all 
2008) later this decade. In addition, characterization 
of several dozens of cold brown dwarf discoveries ex- 
pected from UKIDSS, the Canad a France Hawa i i Tele - 
scope Legacy Survey (CFHTLS; iDelorme et~aTl 120081) , 
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the Wide Field Infrared Explorer ('WISE: llIu~eraII[2008l) 
and the Visible and Infrared Survey Te lescope for As- 
tronomy (VISTA; McPherson et al. 2006) cannot rely on 
space-based spectroscopy. Rather, accurate reproduction 
of ground-based near-infrared spectroscopy, facilitated 
by sensitive instruments such as FIRE, will be necessary 
to study the physical characteristics and demographics 
of these intrinsically faint sources. 

5.3. Is Ross 458C a Planet? 

Finally, we consider whether the properties of 
Ross 458C qualify it as a "planet" (|Basri fe Browrjl2006t 
iSoterl 120061: iMarchil [2007D . In terms of characteristics, 
the planet moniker appears appropriate. Ross 458C is 
gravitationally bound to stellar (hydrogen-fusing) com- 
panions, and itself appears to have a mass below which 
no fusion is expected. Its relative mass is roughly 1:50 
compared to its companions. Its surface temperature is 
well below those of several well-studied, so-called plan- 
ets; e.g., AB Pic b (20001^ K: IBonnefov"etlII[2010h . 
1RXS J160929. 1-210524 b (+800+200 K: lLafreniere eTal 



2010), 2MASS 1207-3932b (1600±100 K; Mohantv et al 
20071 and HR 8799 b (1300-1700 K. IBowler et al.ll2010T) 



Evidence of clouds also lends Ross 458C a decidedly 
planetary-like character. 

However, cosmogony seems to indicate otherwise. 
Ross 458C's 1100 AU projected separation makes it 
highly unlikely to have formed according to the "stan- 
dard" model; i.e., the accumulation of gas and dust 
in a circumstellar ( perha ps in this case circumbinary) 
disk (pollack et al J 119911 . Indeed, if Ross 458C is a 
planet, it has the widest and longest (^50,000 yr) plan- 
etary orbit known, 2-4 times wi der than AB Pic b 
(M » 0.011 Mp, a » 260 AU; IChauvin et all 120051 : 
iBonnefov et~aH 120101) . 1RXS J160929 .1-210524 ( M w 
0.008 Mq, a ~ 330 AU;lLafreniere et al.l l2008ll20Toh . CT 
Cha b (M = 0.016±0.006 AU, a « 440 AU: lSchmidt et al.l 
120011 a nd UScoCTIO 10 6b (M = 0.013±g:Sol> a « 
670 AU; iBeiar et al.ll2008| ). Can planets form at these 
separations? Some disk fragmentation simulations have 
been able to produce very wide (>1000 AU) substel- 
lar companions, albeit with extreme initial conditions 
(e.g., Mdisk/M* » 1: jStamatellos fc Whitworthl 120091: 
iVorobvov fc Basul l2010t see also iBossI 1200611 Planet- 



planet interactions have been shown to scatter sources 
out to orbits 100-1000 AU in size, although they are ex- 
ceedingly rare (<1%) and shor t -lived in simulatio ns (e.g., 
iDodson-Robinson et al.l [20091: iVeras et al.l [2009h . The 
triple nature of Ross 458 may provide an angular mo- 
mentum sink for such a scenario, all owing a planet to be 
scattered out without being ejected ( Holman fc Wiegertl 
[19991 INelsor][2003l: IPierens fc Nelsonll2008fl . All of these 
scenarios require dramatic — and highly fortuitous — 
dynamic interactions to send a planet formed in the re- 
gion of a circumstellar/circumbinary disk to an orbital 
radius 100 times wider than Jupiter's. Of course, a sim- 
pler alternative may be that Ross 458 formed as a "nor- 
mal" hierarchical triple star system, albeit one with very 
low mass components. 

Ross 458C therefore stands as a benchmark not just for 
studies of low-temperature, planetary-like atmospheres, 
but for the definition of the word planet in general. Bet- 
ter characterization of the age of this system and the 



orbital and physical properties of its components will ul- 
timately help us address what parameter space defines a 
planet, and whether "well-defined boundaries, empirical 
verifiabilit y, or logical consiste ncy ...[can] overcome gut 
feelings" (jBasri fc Brownll2006f ). 

6. SUMMARY 

We have obtained FIRE near-infrared spectra of the 
faint companion Ross 458C and the late-type field T 
dwarfs ULAS J1335+1130 and SDSS J1416+1348B. 
Data for the first source confirm its substellar nature and 
identify it as a T8 dwarf, with an inferred log 10 £& z / Lq 
= — 5.62±0.03 and T e /y = 650±25 K based on empirical 
bolometric corrections and evolutionary models. These 
data also show a pronounced K-band peak and relatively 
red J — K spectrophotometric colors, indicative of super- 
solar metallicity and youth that coincide with the prop- 
erties of the Ross 458 system (age = 150-800 Myr, [Fe/H] 
= +0.2 to +0.3). Fits of the data to synthetic spectra 
from SM08 confirm these secondary parameters, albeit 
with discrepant effective temperatures and distances de- 
pending on the presence or absence of condensate opac- 
ity. In fact, cloudy models provide significantly better 
fits to the near-infrared spectrum of Ross 458C and closer 
agreement with its parallax distance, while also providing 
marginally better fits to data for ULAS J1335+1130 and 
SDSS J1416+1348B. The fact that cloud opacity effects 
are more pronounced in the spectra of Ross 458C and 
ULAS J1335+1130 indicates that clouds are particularly 
relevant for shaping the spectra of young, low-gravity 
brown dwarfs, and/or those with supersolar metallicities, 
following trends seen in the warmer L dwarfs. 

Despite the improvement, model fits to brown dwarf 
spectra remain a work in progress. The cloudy mod- 
els "overcorrect" Y- and J-band fluxes and underesti- 
mate absolute fluxes in the case of Ross 458C. Both 
shortcomings suggest that thinner clouds (e.g., f se d — 
3 or 4) may be needed to accurately account for con- 
densate extinction, in line with values inferred for L/T 
trans ition objects (Cushin g et al.l 120081 : iStephens et al.l 
I2009T K Furthermore, the limited sample examined here 
may not be representative of the greater population of 
low-temperature, low surface gravity T dwarfs. A com- 
prehensive modeling effort can more fully explore the 
properties of clouds in these objects, while accurate par- 
allax measurements would provide the necessary empir- 
ical constraints on surface fluxes. Nevertheless, our re- 
sults indicate that by including condensate opacity in 
T dwarf spectral models, more accurate determinations 
of the physical parameters of cold brown dwarfs and 
exoplanets are possible with near-infrared spectroscopy 
alone. 
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